Non-averaged regularized formulations as an alternative to semi-analytical orbit propagation methods

This paper is concerned with the comparison of semi-analytical and non-averaged propagation methods for Earth satellite orbits. We analyse the total integration error for semi-analytical methods and propose a novel decomposition into dynamical, model truncation, short-periodic, and numerical error components. The first three are attributable to distinct approximations required by the method of averaging, which fundamentally limit the attainable accuracy. In contrast, numerical error, the only component present in non-averaged methods, can be significantly mitigated by employing adaptive numerical algorithms and regularized formulations of the equations of motion. We present a collection of non-averaged methods based on the integration of existing regularized formulations of the equations of motion through an adaptive solver. We implemented the collection in the orbit propagation code THALASSA, which we make publicly available, and we compared the non-averaged methods to the semi-analytical method implemented in the orbit propagation tool STELA through numerical tests involving long-term propagations (on the order of decades) of LEO, GTO, and high-altitude HEO orbits. For the test cases considered, regularized non-averaged methods were found to be up to two times slower than semi-analytical for the LEO orbit, to have comparable speed for the GTO, and to be ten times as fast for the HEO (for the same accuracy). We show for the first time that efficient implementations of non-averaged regularized formulations of the equations of motion, and especially of non-singular element methods, are attractive candidates for the long-term study of high-altitude and highly elliptical Earth satellite orbits.Comment: 33 pages, 10 figures, 7 tables. Part of the CMDA Topical Collection on "50 years of Celestial Mechanics and Dynamical Astronomy". Comments and feedback are encourage

Generalization of a method by Mossotti for initial orbit determination

Here we revisit an initial orbit determination method introduced by O. F. Mossotti employing four geocentric sky-plane observations and a linear equation to compute the angular momentum of the observed body. We then extend the method to topocentric observations, yielding a quadratic equation for the angular momentum. The performance of the two versions are compared through numerical tests with synthetic asteroid data using different time intervals between consecutive observations and different astrometric errors. We also show a comparison test with Gauss's method using simulated observations with the expected cadence of the VRO-LSST telescope.Comment: 22 pages, 9 figure

Asymptotic solution for the two-body problem with constant tangencial acceleration

An analytical solution of the two body problem perturbed by a constant tangential acceleration is derived with the aid of perturbation theory. The solution, which is valid for circular and elliptic orbits with generic eccentricity, describes the instantaneous time variation of all orbital elements. A comparison with high-accuracy numerical results shows that the analytical method can be effectively applied to multiple-revolution low-thrust orbit transfer around planets and in interplanetary space with negligible error

Time elements for enhanced performance of the Dromo orbit propagator

We propose two time elements for the orbit propagator named Dromo. One is linear and the other constant with respect to the independent variable, which coincides with the osculating true anomaly in the Keplerian motion. They are defined from a generalized Kepler’s equation written for negative values of the total energy and, unlike the few existing time elements of this kind, are free of singularities. To our knowledge it is the first time that a constant time element is associated with a second-order Sundman time transformation. Numerical tests to assess the performance of the Dromo method equipped with a time element show the remarkable improvement in accuracy for the perturbed bounded motion around the Earth compared to the case in which the physical time is a state variable. Moreover, the method is competitive with and even better than other efficient sets of elements. Finally, we also derive a time element for a null and positive total energy

Accurate analytical approximation of asteroid deflection with constant tangential thrust

We present analytical formulas to estimate the variation of achieved deflection for an Earth-impacting asteroid following a continuous tangential low-thrust deflection strategy. Relatively simple analytical expressions are obtained with the aid of asymptotic theory and the use of Peláez orbital elements set, an approach that is particularly suitable to the asteroid deflection problem and is not limited to small eccentricities. The accuracy of the proposed formulas is evaluated numerically showing negligible error for both early and late deflection campaigns. The results will be of aid in planning future low-thrust asteroid deflection mission

A new set of integrals of motion to propagate the perturbed two-body problem

A formulation of the perturbed two-body problem that relies on a new set of orbital elements is presented. The proposed method represents a generalization of the special perturbation method published by Peláez et al. (Celest Mech Dyn Astron 97(2):131?150,2007) for the case of a perturbing force that is partially or totally derivable from a potential. We accomplish this result by employing a generalized Sundman time transformation in the framework of the projective decomposition, which is a known approach for transforming the two-body problem into a set of linear and regular differential equations of motion. Numerical tests, carried out with examples extensively used in the literature, show the remarkable improvement of the performance of the new method for different kinds of perturbations and eccentricities. In particular, one notable result is that the quadratic dependence of the position error on the time-like argument exhibited by Peláez?s method for near-circular motion under the J2 perturbation is transformed into linear.Moreover, themethod reveals to be competitive with two very popular elementmethods derived from theKustaanheimo-Stiefel and Sperling-Burdet regularizations

Techniques for simulation and control of innovative aerospace systems: numerical models for attitude and orbital dynamics

The PhD thesis deals mainly with special and general perturbation techniques. The main results are the development of a new regularized method for propagating the motion of a particle around an attractive central body and the derivation of accurate analytical formulas to predict the motion in the particular case of constant tangential thrustLa tesi di dottorato tratta principalmente tecniche speciali e generali di perturbazione. I risultati più significativi sono lo sviluppo di un nuovo metodo regolarizzato per la propagazione del moto di una particella intorno ad un corpo centrale di attrazione e la derivazione di formule analitiche accurate per predire il moto nel caso particolare di accelerazione tangenziale costant

EDROMO: An accurate propagator for elliptical orbits in the perturbed two-body problem

EDROMO is a special perturbation method for the propagation of elliptical orbits in the perturbed two-body problem. The state vector consists of a time-element and seven spatial elements, and the independent variable is a generalized eccentric anomaly introduced through a Sundman time transformation. The key role in the derivation of the method is played by an intermediate reference frame which enjoys the property of remaining fixed in space as long as perturbations are absent. Three elements of EDROMO characterize the dynamics in the orbital frame and its orientation with respect to the intermediate frame, and the Euler parameters associated to the intermediate frame represent the other four spatial elements. The performance of EDromo has been analyzed by considering some typical problems in astrodynamics. In almost all our tests the method is the best among other popular formulations based on elements